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ABSTRACT 

We present for the first time predictions for UV line emission of intergalactic and cir- 
cumgalactic gas from Adaptive Mesh Resolution (AMR) Large Scale Structure (LSS) 
simulations at redshifts 0.3<z<1.2, with a specific emphasis on its observability with 
current and near-future UV instrumentation. In the three UV transitions of interest 
(Lya, OVI and CIV) there is a clear bimodality in the type of emitting objects : 
the overwhelming majority of the flux stems from discrete, compact sources, while a 
much larger fraction of the volume is filled by more tenuous gas. We characterise both 
object types with regard to their number densities, physical sizes and shapes, bright- 
nesses and luminosities, velocity structures, masses, temperatures, ionisation states, 
and metal content. Degrading our AMR grids to characteristic resolutions offered by 
available (such as FIREBall) or foreseeable instrumentation, allows us to assess which 
inferences can be drawn from currently possible observations, and to set foundations 
to prepare observing strategies for future missions. In general, the faint emission of 
the IGM and filamentary structure remains beyond the capabilities of instruments 
with only short duration exposure potential (i.e. stratospheric balloons), even for the 
most optimistic assumption for Lya, while the yet fainter metal line transitions (OVI 
and CIV) for these structures will actually remain challenging for long duration expo- 
sures (i.e. space-based telescopes), mostly due to their low metallicities pushing them 
more than three orders of magnitudes in brightness below the Lya radiation. For the 
bright, circum-galactic medium (CGM) the situation is much more promising, and it 
is foreseeable that in the near future we will not only just dectect such sources, but the 
combination of all three lines in addition to velocity information will yield valuable 
insight into the physical processes at hand, illuminating (and discriminating between) 
important mechanisms during the formation of galaxies and their backreaction onto 
the IGM from whence they formed. 

Key words: Proper keywords here. 
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Both inventor ies of the cosmic bar y ons budget 

l|Fukugita et al.l Il998l : iFukugita fc Peebles! |2004l). and 
cosmological simulations ( Cen fc Ostrikej 199E ; Dave et"aLl 
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1 19991 . I2OOII : ICen fc Fanj I2OO6I ) which predict the de- 
velopment of a warm-hot phase of the intergalactic 
medium (IGM) at low redshift have motivated an in- 
tensive search for baryons i n a p h ysical state not we ll 
explored yet (e 
iDanforth fc Shul 



Lehner et all (I2OO1I): 



or baryons m a p r 
g. iBregmanI ll2007r): 

l1(l2008ll :lTripp et all ., . 

l|2008l '): IDanforth et al.l (|201Cl ) and references therein) 



l (|2008h : lThom fc ChenI 



Beyond the baryon accounting issue itself, this approach 
is providing new insights into important aspects of cos- 
mic evolution, specifically about the low fraction of baryons 
(about 6 %) that have made their way into stellar popula- 
tions of galaxies, the interdependence between the IGM and 
galaxy evolution, and the crucial role of baryons circulating 
in and out galaxies. 

The exploration of matter in the IGM at densities 
below the mean density is naturally performed through 
quasar (or other bright background sources) absorption-line 
studies, and has a long history of spectacular successes. 
At low redshift, in the domain of tracers of the Warm-Hot 
intergalactic medium (WHIM), this is now culminating 
in observations with the Cosm ic Origins Spect r ograp h 
(COS) on board of HST (e.g. iNaravanan et al.1 (|201lf ): 
ISavage et all l|201ll ). 

Absorption line studies are, however, traditionally 
limited by a number of factors, such as the need for 
bright background sources, the difficulty of exploring 
perpendicular to a line of sight, and the difficulty of finding 
background sources for investigating a specific object. In 
addition any improvement in sensitivity requires an increase 
of c ollecting area. Emiss io n line observ ations, as discussed 
by iHogan fc Wevmaiml (|l987l ) and ICould fc Weinberg! 
l|l996tl , have therefore been considered as an alternative 
approach. They are not affected by the issue of background 
sources and, depending on the observed solid angle, are 
less demanding than absorption techniques i n ter ms of 
collecting area. In this vein, e.g. iRauch et ahi (120081 ) have 
performed a very long exposure with VLT in order to 
detect Lyman a emission from matter at the transition 
from ionized gas to neutral, self-shielded gas. 

Exploiting a variety of simulations, predictions have 
been made for the possible signatures of structure formation 
and 'feedback processes' in terms of emission lines over a 
range of redshifts. Focusi ng on hydrogen a s the n a turall y 
most abundant spec ies, iFurlanetto et al.l (|2004l . |2005| ): 
iKoUmeier et al.l (|2010l ) and references therein, predict Lya 
emission from IGM gas at ^ 2 ^ 5, and discuss various 
mechanisms for its origin (reradiation of absorbed ionising 
photons from the metagalactic UV background; radiative 
cooling of shocks during the assembly of gravitationally 
bound objects; interaction of strong winds with the sur- 
rounding IGM). Furthermore, these authors also stress 
the difficulty to discriminate between these, in addition 
to the compli c ation that galaxies them s elves emit in Lya. 
iBertone et al] (|2010D and lBertone et al.l (|2010l ') trace a large 
suite of metals in different ionisation stages through their 
SPH simulation in combination with state of the art gas 
cooling prescriptions to track their emission both in the 
soft-X rays as well as in the (restframe) UV over a similar 
redshift range as above(0<z<5). Of particular interest 



are their notions that the transitions in CHI, CIV, SilV 
and OVI are bright enough (>1000 photons/(s cm^ sr) to 
allow for a detection with current technology ('in the near 
future') in gas of moderate overdensities [5 ^ 100) that has 
been enriched to metallicities Z> O.lZsoiar- They caution, 
however, that hence these emission lines are biased tracers 
of the baryons, but still provide good tools to detect gas 
cooling onto or flowing out of galaxies. 

In this paper, we follow a similar line of investigation 
for Lya, CIV and OVI emission over the redshift inter- 
val 0.3<z<1.2, but go one step further in the detailed 
effects of observability with current and/or near-future 
instrumentation in the UV. We clearly highlight the 
limitations posed by the detectors' spatial resolution, and 
stress the importance of taking the velocity structure 
of the emitting gas into account. We run source finding 
software in simulated data-cubes on grids resembling a 
typical instrumentational set-up, and are for the first time 
able to characterise potentially detectable structures with 
respect to their brightness, luminosities, masses, surface 
and volume densities, sizes and morphologies. There is a 
clear delineation (in Lya) between the very diffuse and 
extended filamentary structures tracing the true IGM, and 
'haloes' surrounding the bright and massive nodes of the 
Cosmic Web, which are more compact (< 200 kpc/h in 
diameter) and much brighter than the filamentary sheets. 

The paper is organised as follows. We briefly discuss the 
simulations in Section 2, highlighting the methods to calcu- 
late the gas emissivities in l2.2l and add a description how to 
arrive at an observable data cube in l2.5l In section Owe take 
a first look at these cubes in all three emission lines we are 
primarily interested in (Lya, CIV, and OVI), and discuss 
general features of the emission in terms of brightness dis- 
tribution functions and spatial extent, while sections [5] and 
2] provide a closer look at the extended, filamentary struc- 
tures and the compact, bright sources, respectively. Before 
concluding in section [71 we discuss the implications of our 
findings in terms of observing strategies and optimal tech- 
nical instrument specifications in sectional 
Throughout the paper we use a ACDM cosmology with the 
results of WMAP5 (see next section for detailed parame- 
ters). 



2 SIMULATIONS 

2.1 Cosmological simulations 

The high-resolution cosmological simulation investigated in 
this article is part of the BINGO! simulatiorQ suite. It 
was performed at the CINES computational center with 
the RAMSES Adaptive Mesh Refinement (AMR) code 



^ BINGO! is a 4-year project of the "Agence Nationale de la 
Recherclie" (ANR), started in 2009. The consortium of three part- 
ner institutes (Marseille, Lyon, and Paris) specifically investigates 
both theoretically and observationally the physical processes of 
structure and galaxy formation with an emphasis on their inter- 
actions with the IGM and CGM. 
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ijTevssierl l|2002h V Initial cond i tions were set using the MP- 
GRAFIC code (iPrunet et al.1 hOO^). a parallel implemen- 
tation of GRAFIC (jBertsch ingei^ ('2001I)). and made use of 
the IfORIZON project white nois<jf]. They assume the fol- 
lowing cosmological parameters, consistent with the results 
of WMAP5: n^(0) = 0.26,f7A(0) = 0.74, /i = 0.719, D.b = 
0.044, as = 0.796, n, = 0.96. 

The comoving box length is 100 Mpc/h , the number 
of dark matter particles is 512"^, the number of cells at the 
coarse level is 512'' as well, and the number of refinement 
levels is 7. As a consequence, the dark matter particle mass 
is TUp = 4.42 X 10* Mq/Zi, and the spatial resolution at the 
highest refinement level is roughly constant and equal to 
1.53 physical kpc/h. 

RAMSES computes the evolution of dark matter parti- 
cles, gas cells, and star particles within the AMR hierarchy 
of grids. Particles are evolved using a PM method with a 
multigrid Poisson Solver. Euler equations for the gas are 
solved with a second-order unsplit Godunov method. We 
adopt the UV background from Haardt and Madau (1996) 
with a zero percent escape fraction from galaxies below the 
Lyman limit. For instance, the flux at z=0.75 and at 912 A 
is 11.75 photons cm~^s~^sr~^ A~^. The net cooling rate in 
the simulation is dependent on the density and metallicity 
of the gas, using a parametrization based on cooling curves 
calculated with Cloudy for a gas in ionization equilibrium 
with the adopted UV background. 

It is worth noting that most of the work about UV 
emission lines of the WHIM/CGM was done with SPH 
ijFurlanetto et al.l (|2003) : iBertone et al.l (|2010l )l or fixed 
grid codes (jCen fc Fand (|2006')). whereas the modelisation 
of mildly-overdense shock-heated regions might be subject 
to caution with these Lagrangian or low-resolution codes. 
While Smith et al. (2011) also use an AMR code (ENZO), 
they are focusing mostly on characterising the WHIM and 
OVI absorption rather than emission from CGM and IGM 
gas as in our case. 

Star formation and feedback follow the methodology de- 
scribed in Dubois and Teyssier(2008, 2010). Star formation 
occurs in cells denser than a density threshold of nu = 0.1 
at/cm^, at a rate proportional to the gas density with an 
efficiency of 2% per free-fall time. The effect of type II SN 
explosions is taken into account by releasing within a ra- 
dius rsn = 15 kpc from each stellar particle, lOMyr after its 
formation, (i) 10% of the mass formed, (ii) lO^lerg per SN 
event, half in the form of kinetic energy and half in the form 
of thermal energy, and (iii) metals, assuming a global yield 
of 0.1. Metals are then advected as a passive scalar enriching 
the ISM, but also the CGM and WHIM of interest here in 
this article. Finally, a polytropic equation of state of index 
7 = 5/3 is used to provide pressure support to the dense gas 
(tih > 0.1 at/cm^) and avoid artificial fragmentation. 



2.2 Gas emission mechanisms accounted for by 

the simulation (photoionized and collisionally 
excited gas) 

The estimate of gas line emissivities follows naturally 
from the simulations through the use of grids indexed by 



hydrogen density and temperature, assuming the same 
uniform ionizing background as adopted for the simulations 
and hence being self-consistent regarding this detail. We 
have used the c8.00 version of the spectral synthesis 
code Cloudjlfl, last described by Ferland et al. (1998), to 
compute the line emissivities of a plasma submitted to 
an incident continuum. The line emissivities of interest in 
this paper are Lya at 1216 A, OVI at 1032/1038 Aand 
CIV at 1548/1551 A. Grids have been generated at the 
redshifts z=0.35, 0.76 and 1, and for each redshift, the 
incident continuum used in Cloudy is the same ultra-violet 
background used in the simulation at the same redshift. 
Practically, these specific redshift values have been adopted 
in order to place the redshifted lines in the unique narrow 
spectral UV window at 2000 A accessible from stratospheric 
balloon-borne experiments such as FIREBall; in the fol- 
lowing we will explore the redshift influence since our goal 
is to understand gas emission over the redshift range 0.3 - 1 . 

Grid points are computed for temperatures between 
3 X 10^ K and 3 x 10* K every Alog{T/K) = 0.01 and 
for densities between uh = 1.778 x 10~* and uh = 1 
cm"'' every Alog{nH /cm~'^) = 0.5. Two sets of grids 
are calculated, one assuming the plasma to be in pho- 
toionization equilibrium with the UV background and 
another one assuming coUisional ionization equilibrium. 
The grids are computed at solar metallicity and the emis- 
sivities of each gas cells is rescaled according to their gas 
metallicity. We use the "no induced proc esses" of Cloudy 
for a n explanation for this choice see e.g. iFurlanetto et al] 
2004h and the following treatment of the photon pumping) . 



2.3 Self-shielding 

Our simulation assumes a uniform ionizing background 
which heats the gas as if it were optically thin everywher^fl. 
This introduces two potentially important limitations to 
our predictions of gas emissivities: (i) the ionization state 
of self-shielded gas is not accurately described (the neutral 
fraction is under-estimated), and (ii) the temperature of 
self-shielded gas may be over-estimated (heating from the 
UV background is active everywhere). 

Self-shielding of gas to cosmological (inhomogeneous) 
ionizing background is an extremely complex question 
which can really be addressed only with simulations which 
jointly solve for radiative transfer and hydrodynamics. 
There has been substantial work to accomodate for that 
and describe a posteriori the effect of self-shielding, ei- 
ther in the perspect ive of absorbers (e.g. iS_chavc 2001a^; 
IPont zen et aP l2008l. and r eferen c es therein), "cooling" 



radiation l|Furlanetto et al. 12004 Faucher-Giguere et al.l 
|2010| . e.g.), fluorescence ( Kollmeier et al.r2010l '). or of 



www.projet-horizon.fr 



course in the perspective of high-redshift reionization (e.g. 



http://www.nubIado.org 
^ In practice, regions with njj > O.lcm"^ (and T < lO'^K), are 
described with a polytropic equation of state, and are thus effec- 
tively seff-shielded since their temperature is solely a function of 
their density. 
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Figure 1. The relative flux distribution for the simulated 100 Mpc/h box in the log T - log n^f bins (coloured phase-space diagram), 
and the Ly a volume emissivity (black solid lines, in units of erg/(sxcm'^)) as function of hydrogen density and gas temperature using 
only photo-ionisation equilibrium (PIE, left panel) or coUisional ionisation equilibrium (CIE, right panel) assumptions at a redshift of 
z=0.75. In addition, the location of the cuts used to delineate various scenarios for the self-shielding treatment is indicated by the dashed 
lines plus the two rectangular boxes on the lower right (for details, see text). It is clear that for Lya emission the exact location of those 
cuts plays a vital role in determining the total emission, as gas with high density and T~ 10* K dominates the emission. Note that in 
this and the following plots, the choices of redshifts stem from our restriction to a constant wavelength coverage over a well-defined UV 
regime of a putative instrument, in order to be realistically illustrative. 



lAubert fc Tevssierll2010l . and references therein) . Such mod- 
eling i s beyond the scope of th e present paper. Instead, fol- 
lowing |^rlanettoIe£^l] (|2004l ). we chose to bracket the real 
solution with a series of simple assumptions : (i) a crite- 
rion to decide which gas is self-shielding (depending on the 
emission line of interest), and (ii) crude estimates of how 
self-shielding gas should emit. 

Self-shielding thresholds 

Because we are interested in emission of IGM or CGM 
gas only, we first discard the gas with nn > nisM = 0.1 
atoms/cm'^. In our simulations, this gas is s tar-forming ISM. 
In the approach bv lFurlanetto et al.l (|2004l ) e.g. it is decided 
that gas with nn > nss and T < Tss is self-shielding. Here, 
Tss is a temperature above which the emitting species (e.g. 
Hi) is mostly ionized by collisions (e.g. at Tss ~ 10*'^, H is 
mostly in the form of Hll, so that the gas becomes transpar- 
ent to Hi-ionizing photons). Above Tss, the emitting species 
cannot self-shield, and photo-ionization becomes dominant, 
nss is some typical density at which self-shielding should be 
effective. In reality, self-shielding is a column density effect. 



and using a cut on density only is bound to be a poor proxy. 
Here, we have in addition chosen another, more physically 
motivated approach. We assume that only gas at a suffi- 
ciently high pressure self-shields, and take the thermal pres- 
sure (P/k ~ nT) as a proxy for the total pressure that guar- 
antees hydrostatic equilibrium (including e.g. magnetic pres- 
sure, cosmic ray pressure etc.). Empirically, it is found that 
a threshold of P/k = 155 cm^ K r esults in an HI density of 
p = 6.1 X IQ^hMQMpc^ (at z=0, [Popping et"al] l|2009l )). a 
value that is similar to the value observed in Zwaan et al. 
(2003) for their self-shielding thresholds. Furthermore, only 
gas whose radiative recombination time (rrec ~ x rie) 
is shorter than the sound crossing time (ts = R/Cs, with 
the sound speed Cs ~ T"'^ and the relevant length scale 
R) can recombine, and hence self-shield. A priori, it is not 
clear which length scale should be chosen, here we have as- 
sumed two values (R=1.0 and R=0.1 kpc). The densities 
and temperatures delineated by these boundaries are within 
the ranges given in th e literature (e.g . Zwaan et al. (l2003h: 
IWeinberg etai] (|l997l ): IWolfire et aLr(|2003l ) and iPelupessvl 
l|2005M . For the metal line transitions similar arguments 
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Figure 2. Same structure as in Figure [T] now sliowing OVI at a redsliift of z=1.0. Note how liere (and in the following plot for CIV) 
the area of main emissivity shifts to much higher temperatures (above 10^ K), and thus the exact location of the self-shielding cut does 
not play an important role anymore for determining the total flux. 



hold. Further details of this me thod to deal w i th th e self- 
shielding limits are listed e.g. in iPopping et al] l|2009l ). 

These thresholds are summarised in Table [T] and graph- 
ically represented by the dashed and solid lines in Figs. 1-6. 
Note that while we have analysed the simulation output sub- 
jecting it to a larger number of such cuts (which themselves 
are not part of the modelling, yet simply a-posteriori im- 
posed cuts on certain grid cells, whose emission we either 
modify or ignore completely), we will focus in the following 
on the two most extreme cases for bracketing (subsequently 
called CUT 7 and CUT 8), in order to simplify the presen- 
tation. 



Emissivity of Self-shielding gas 

We use two simple scenarios to compute the emissivity of 
the self-shielding gas. Either we assume this gas to be in 
coUisional ionization equilibrium (CIE) at the temperature 
and density provided by the simulation and use CIE tables 
to predict its emissivity. Or we argue that had the simulation 
taken self-shielding into account, the gas would be at much 
lower temperatures than predicted, and would simply not 
emit even in CIE. In that case, we set the emissivity to zero. 



2.4 Other emission processes: photon pumping for 
Lya radiation 

There is deliberately no light emission associated with 
the formation of stars in this simulation. Howev er, the 
UV background used here (|Haardt &: Madaul 1200 ll ) , does 
include on top of the QSO emission an integrated contri- 
bution from stars, but only above the Lyman limit (ie. 
zeroed escape fraction below that limit). Due to the large 
cross-section for absorption of hydrogen gas at Lyman line 
wavelengths, continuum photons from the non-ionizing 
background can be scattered when they are redshifted to 
the Lyman resonance lines. This process, often coined 'pho- 
ton pumping', has not been accounted for by our simulation. 

With an isotropic background, as assumed here, this 
pumping does not have a net effect for Lya emission, as it 
simply consists of a redirection photons in a medium that 
is permeated by an already isotropic radiation field - the 
only true additional contribution may arise from the con- 
version of hi gher orde r Lym an line photons into Lya pho- 
tons. iFurlanetto et al] (|2004i ) estimate this contribution to 
remain below the level of a few percent. Hence, our calcu- 
lations of the gas emissivity explicitly avoid the effects of 
photon pump ing by turning th i s featu re off in CLOUDY, as 
suggested bv IFurlanetto et al.l (|2004l ) . Note, however, that 
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Figure 3. Same structure as in Figure [T] now sliowing CIV at a redshift of z=0.35. 



the situation changes locally, if the UV radiation field is 
anisotropic and potentially of different spectral shape than 
the general background - as it can be imagined near UV 
bright galaxies. Properly modelling the photon pumping 
contribution (or actually tracing the continuum photons re- 
sponsible for it on their paths through the simulated cubes) 
is beyond the scope of this paper, but we have performed 
calculations with simplified assumptions to judge the impor- 
tance of it (cf. 6.2). From these, we conclude that indeed en- 
hancements of the fiux over the emissivities calculated here 
of an order of magnitude are possible in the proximity of 
UV bright sources. 



2.5 Simulated Observations 

We construct mock data cubes from our simulations in two 
steps: first, we build a pre- observation light-cone at the full 
AMR resolution, and second, we convolve this light-cone 
with instrumental effects. 

Pre-observation light-cones 

The typical geometry of the light-cones we wish to produce 
is that of a pencil beam, i.e. a cone with a radial extent much 
longer than the size of our simulation (typically 5 times) 
and an angular extent probing scales less than a tenth the 
size of our simulation at most. Another point in considera- 
tion is that we are interested in emission from diffuse inter- 



galactic gas, which can produce extended signal on scales up 
to a few tens of Mpc. It is thus important to preserve the con- 
tinuity of the gas densi ty field in the lig ht-cones. We use the 
MOMAF software (jBlaizot et al.ll2005l ) to cut a hght-cone 
out of our simulation. MOMAF takes care of replicating our 
(periodic) simulation at will and to cut light-cones in any 
direction relative to the axis of the box. In order to preserve 
continuity, we chose not to use the random tiling technique 
and to build light-cones out of one simulation output only. 
This avoids both re dshift and replication discontinuities (see 
iBlaizot et al.|[2005l . for details). We then chose lines of sight 
by carefully selecting angles for light-rays in order to mini- 
mize replication effects. 

In practice, we build three different types of data cubes 
(see below, one for each emission line of interest) For this, 
we use the snapshots, i.e outputs at z ~ 0.35 (resp. 0.75, 
1.04) for Civ (resp. Hi, Ovi). We have checked that our 
conclusions are not changed if we use slightly different snap- 



® Here, and in the following, we use the term 'data cube' for a 
specific output of the simulation, where we have taken a part of 
the whole simulated volume and processed it in different ways. 
For some of the latter analyses, e.g. we have created regularly- 
gridded data cubes, whereas other aspects of the analysis may 
rely upon the original AMR gridding. As such, the terms 'cone' 
and 'cube' are for practical reasons equivalent in their meaning. 
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Table 1. Different models for the emissivity of the gas in the light cones. Generally speaking, we divide the parameter space in loci 
that exhibit either strictly no emission, or an emissivity is calculated by applying CLOUDY with either photo-ionisation equilibrium 
(PIE) or coUisonal-ionisation equilibrium (CIE). The number indicat ed for each model corresponds to the one used in the different plots. 
Models 2 and 3 follow closely the procedure and values suggested bv Furlanetto et al.l {2004), while models 4 and 5 allow for a ten times 
higher limiting density threshol d for CIE or PIE. Mo dels 6-9 involve the pressure and timescales for recombination and sound crossing, 
as described in the text and in [Popping et al. I l l2009t) , for different size scales (R=0.1 and R=1.0 kpc). The units for the density are 
[n]=cm~'^, for the temperature [T]=K, and [P/k]=cm~'^K. As figures 1-6 show, the cuts labeled 7 and 8 are the most severe and lenient, 
respectively, for the line emission, and hence we use them to bracket the intermediate cases. 





Cut 




Zero emissivity 




PIE emissivity 


CIE emissivity 


SF cut 


1 




nu > 0.1 




nn < 0.1 




-l-CIE 


2 




nn > 0.1 




log T>4.5 OR 


log T< 4.5 AND 








nn < 5.1 X 10-3 (Lya) i 




nn ^ 5.1 x 10-3 (Lya) 




-I- cut 


3 




nu ^ 0.1 OR 




log T>4.5 OR 








(log 


: T< 4.5 AND nn ^ 5.1 x 10" 


3(Lya)) 


nn < 5.1 x 10-3 (Lya) 




-l-CIE xlO 


4 




nn ^ 0.1 




log T>4.5 OR 


log T< 4.5 AND 








nn < 5.1 X 10-2 (Lya) 2 




nn Ss 5.1 X 10-2(Lya) 




+cut xlO 


5 




nn > 0.1 OR 




log T>4.5 OR 
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"(Lya)) 


nn < 5.1 X 10-2(Lya) 




R=l kpc -1- CIE 


6 




riH > 0.1 




P/k <258 OR Tree > Ts ^ 
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P/k >258 AND Tree < Ts 
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nn > 0.1 




P/k <258 OR Tree > Ts 


P/k >258 AND Tree < Ts 


R=0.1 kpc -l-cut 
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P/k >258 AND Tree < Ts 
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1 The values change with different redshifts, hence : nn > 7.1 X 10-3(0^/) and nn > 2.7 X 10-3(C/V). 
2 n.H > 7.1 X 10-2(Oy/) and nn > 2.7 x \Q~'^{CIV). 
3 For OVI at Z'-l.O, we apply P/k(z) = 356, and for CIV at Z'-0.35 P/k = 135 as dividing line. 



shots, or if we use them all and break the continuity of the 
density field. 

In the outputs, each gas element (cell) has a luminosity 
distance, an apparent redshift (taking into account the gas' 
peculiar velocity), and angular coordinates. 

Observed data cubes 

As an instructive example, we chose to mimick the in- 
strumental characte ristics of FIREBall (|Tuttle et al.ll201ol : 
iMilliard et "alll2010l ). We will discuss alternative choices of 
angular or spectral resolutions in Sec. 5, but note here 
that current or near-future UV missions are of quite sim- 
ilar layout. While focusing on the FIREBall spectral cov- 
erage forces certain specific redshift choices for the three 
lines in question, it does also allow us to check on differ- 
ences regarding the cosmic evolution over that redshift span 

(0.35 < «em < 1.1). 

In order to produce observed data cubes from of our pre- 
observation catalogs, we proceed in 3 steps. First, at the full 
AMR resolution, we compute the luminosity of each cell ac- 
cording to the various models described in Sec. 12.21 Second, 
we derive each cell's fiux contribution to a regular grid (in 
Q, (5, and z) with better resolution than the typical angular 
and spectral point-spread function. Third, we then convolve 
these in the spectral as well as the spatial dimensions to 
arrive at the desired final resolution (in angular direction 
we apply a sharp disk of an 8 arcsec diameter, and in the 
spectral direction we convolve with a Gaussian of FWHM = 
2.25 pixels ~ 0.56 A), and finally save the convolutions on 
a regularly gridded data-cube. 

We have generated 300 of such fields of view with 
an approximate 3x3 arcmin^ FOV (similar to FIREBall's 
FOV) in order to assess the effects of cosmic variance 
(e.g. for the construction of the PDFs). Furthermore, 
for instructional purposes, we have specifically chosen 3 



sightlines (for each transition and redshift) with a larger 
FOV : one towards a field containing a known very massive 
Dark Matter structure to represent a best case-scenario 
for 'bright' sources, one which contains an example for a 
bright, large filamentary structure, and another one as the 
'best' compromise between those two criteria. These fields 
cover a field-of-view (FOV) of fixed size (900 x900 arcsec^), 
and fixed spectral coverage (270 A), thus creating cubes of 
different volumes depending on the specific redshift. The 
selection of these sightlines, in fact, used PDFs (see below) 
from fields of smaller FOV, in order to determine their 
brightness distributions. 

Figures 1 to 3 show the phase space diagrams of the 
simulation, each cell weighted by its luminosity in the Lya, 
OVI or CIV transitions. Note how dominant the high den- 
sity areas (log n^f ^ —2) are in these, as emissivity (denoted 
by the solid, curved lines) scales with n|f . It is immediately 
obvious that the treatment of self-shielding is very impor- 
tant for the Lyof estimates, while it plays less of a role for 
the metals, which emit most of their light from gas at tem- 
peratures ^ 10^ K. 



3 GENERAL FEATURES OF THE LINE 
EMISSION 

Even a casual visual inspection of the line emission inside 
the simulated cubes (cf Fig. [4]) reveals a clear bimodality 
in the types of emitting regions. There are on the one 
hand bright sources, spatially confined to small areas, but 
extended in wavelength due to the gas being spread out 
in velocity space. In the vicinity of these, and connecting 
such regions, there are filamentary structures, which are 
much fainter, more extended spatially and affected less 
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by the velocity dispersion of the gas. We will discuss the 
specific properties of these two different classes of emitters 
in the next two sections, but focus here first on general 
implications. 

Figures (5] [6l and [7| show the probability distribution 
functions (PDF) of the voxels' emission in Lya, OVI and 
CIV, weighted by volume. The solid lines represent the 
PDFs based upon the full AMR resolution, while the 
dashed lines are PDFs calculated after resampling of the 
data cubes with a grid of 3 arcseconds spatially and 0.25 A 
spectrally, subdividing the simulated cube smoothed with 
the resolution of 8 x 8 arcsec^ and 0.56 A as described 
earlier (resembling the FIREBall instrument specifications) . 
Note that unlike e.g. the PDFs pre sented bv lBertone et al] 
\2Ql &i and ' Furlanetto et al.l (j2004h these are PDFs based 
upon the full 3D voxel distribution, i.e. we do not collapse 
along one axis and derive the distribution on the resulting 
2D images of such slices. Hence, a comparison of the curves 
is not straightforward. The various cuts to deal with the 
self-shielding gas are indicated by the different colours. 
Obviously, the two extreme distributions are achieved when 
the cuts are most lenient or most severe in excluding gas 
located in the parameter space that exhibits the highest 
emissivity. The former happens when we choose to apply 
the CIE case for R=0.1 kpc (denoted as example 8 in Table 
1), the latter is achieved by neglecting all emission for 
R=1.0 kpc (example 7). In the following, and in various 
graphs, we will refer to these two extreme scenarios as 
CUT7 and CUTS, respectively. 

A decrease in resolution has a dramatic effect on the 
brightest areas, but leaves the distribution unchanged for 
surface brightnesses less than ~ 1 Cljf|. The upward shift 
of the plateau ranging from about 2.0 ^ log CU ^5 7.0 at 
the AMR level to a similar plateau with a lower cutoff 
high surface brightness, but a higher volume fraction at 
the lower resolution clearly indicates that the brightest 
pixels are strongly 'clustered', i.e. form collectively the 
bright, spatially compact sources. These bright regions 
are rare : even for the resampled cubes they do not cover 
more than 0.1% of the volume, but carry almost all of the 
light emitted in all three lines. The various cuts dealing 
with the self-shielded gas allow us to bracket the expected 
emission. Note that in the case of Lya the cut-off values 
for the brightest sources change by a factor of 100, a 
direct reflection of the fact that this is the location of the 
parameter space where the emissivity for Lya peaks (cf. 
Fig. [ij. In contrast, as expected, the exact treatment of the 
self-shielded gas has very little effect on the bright end of 
the surface brightness distributions for the OVI and CIV 
line emission. 

Figure |8] shows two PDFs for Lya at the redshifts of 
z=0.75 and z=0.35. Note that there is by no means a simple 
(l-l-z)* evolution at all, as may be expected from sheer cos- 
mological expansion effects (meaning that the lower redshift 
values should exceed the higher ones by a factor of about 
3), but a more complex evolution. The faint emission peak 

6 1 CU = 1 photon s-lcm-2sr-lA-l 



e.g. becomes even less bright, mostly due to the change in 
the ionising radiation field weakening towards lower z. The 
brighter regions beyond the knee in the distribution, how- 
ever, do become brighter at the lower redshift. Hence, the 
redshift evolution of the surface brightness distribution is a 
complicated function, whereby geometric cosmic expansion 
and true evolutionary source effects combine in a non-trivial 
fashion. In any case, the evolution expected from redshift 
changes surveyed here is probably small compared to the 
uncertainties introduced by the self-shielding treatments. 



4 BRIGHT, COMPACT SOURCES/CGM 

4.1 Introduction to the source finder 

The simulations with the reconstructed distribution of 
Lya, CIV and OVI contain bright objects with extended 
environments connected by filaments spanning up to the 
whole spatial scale of the simulation output. Features in 
the simulated outputs only represent circum- and inter 
galactic g all the star forming gas and gas in the 

galaxies has been removed using thresholds as described 
in 2.2. To extract objects from the simulated c ubes, we 
have used the source finding algorithm Duchamp (jWhitind 
l2008h . an algorithm that has specifically been developed 
to detect objects in 3D radio observations that typically 
have dimensions of RA, DEC and velocity. The software is 
however very fiexible and can be used for any kind of three 
dimensional datasetQ 

The simulation outputs tailored to the spectral and 
angular resolution specifications of FIREBall (see earlier) 
were used as input for Duchamp. Within Duchamp, a 
parameter input file constrains the properties of a possible 
detection. For all objects a requirement has been set, 
that detections consist of at least two adjacent voxels 
in the spectral direction. Because of the oversampling in 
the spectral direction by about a factor two, this is the 
minimum size an object can have. 

Within Duchamp objects were sought using two differ- 
ent thresholds. The first threshold represents the minimum 
peak fiux a feature should have to be accepted as a compact 
object. In the reconstructed Ly-a cubes compact objects 
can be distinguished by their large peak fiux. The brightest 
regions in the filaments have typical values of the order of 
~ 10 CU, therefore the first threshold for Ly-a is chosen at 
100 CU. Once the location of bright and compact features 
has been determined, the objects are grown in size by 
adding adjacent voxels until a second (lower) threshold 
has been reached. By using two different thresholds we 
can isolate compact objects down to a relatively low 
flux threshold which is below the peak flux of extended 
filamentary structures. For Ly-a the second threshold is 
chosen at 10 CU. Although detections could be extended 
to lower flux thresholds, this would increase dramatically 
their size and include parts of filaments or possibly even 



The algorithm is decribed and can be obtained here : 
|http://www.atnf.csiro.au/peopIe/Matthew.Whlting/Duchamp/| 
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Figure 4. Examples of the simulation output for Lya line emission at z=0.75. The lower right inset shows a part of an AMR lightcone 
with an opening angle of 15x15 arcmin^ in the x-y direction and a length in the z direction of ~100 A corresponding to a structure in 
(physical units) of 6.6 X 6.6 X 300 Mpc. The web-like structure with long, elongated filamentary bridges connecting the bright nodes (red 
labelling) is clearly distinguishable. The top panel shows the same volume, but now both degraded to a resolution resembling currently 
available instruments (i.e. about 8 arcseconds in spatial, and 0.5 A in wavelength direction) and with peculiar velocities added to each 
cell. Note how these tend to elongate the bright sources in z-directions, while leaving the filamentary structures almost untouched. Also 
note how the bright sources are clearly clustered, and not randomly distributed. 



merge them, without adding significant amounts of flux. 



For each transition the threshold fluxes were deter- 
mined empirically by carefully inspecting the cubes, the 
exact values are given in table [21 The sharp edges of the 
source regions make the number of the detected sources 
largely insensitive to the threshold values, as it is possible 
to choose a value for the second threshold that allows for 
picking up a large fraction of the total light inside the 
simulated volume while ensuring that there is little source 
merging. The same flux thresholds have been used for 
all self-shielding cuts, since the selected cuts only affect 
the densest regions, and even the most severe such cut 
does not prevent detection with the parameters chosen here. 



Table 2. Flux thresholds used within Duchamp to find detections 
in the simulated cubes of Ly-«, CIV and OVI. 



Element 


threshold 1 [CU] 


threshold 2 [CU] 


Ly-« 


100 


10 


OVI 


1 


0.5 


CIV 


0.5 


0.05 



4.2 Source properties : Space density, 
luminosities, sizes, shapes, spectral 
information 

The number, and hence space density (not the flux) of those 
sources in Lya is very robust and almost independenl|f| of 
the self-shielding cut criteria and the specific parameters we 



" For the harshest cut criteria, a select few (in the extreme case 
4 out of > 1000) Lya objects fall below the minimum detection 
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Figure 5. The probability distribution function (PDF) of the voxels' line emission in Lya for z~0.75, weighed by volumc(left panel; i.e. 
we simply count all the cells' volume per certain brightness bin and divide by the total survey volume in order to derive the fraction per 
bin). Note that this is hence a 'tr ue' pdf in the sense tha t w e are using the compl ete 3D cube, and do not collapse along one dimension 
(in contrast to e.g. the PDFs by iFurlanetto et al.l ||2004| ) or iBertone et al] 1I2OIOI) . whose results are based on the pixels of an image). 
The sudden change in slope when reaching the plateau around a few CUs from the low surface brightness side is indicative of reaching 
a different population, which can - as detailed in the text and Figs . [TTI and [121 - be equated with bright compact sources. The effect of 
resampling from the AMR refinement level (solid lines) to the observationally dictated resolution (dotted lines) is a dramatic loss of the 
highest surface brightness objects. Note also that the choice for different treatment of self-shielding results in factors of 10-100 for the 
brightest spots. Bright sources are rare, but carry almost all of the light : 10~* of the volume after resampling. The most important 
aspect regarding the resampling, however, is that the very brightest spots still remain within the range of capabilities of an instrument 
with current technology, as discussed in section [6] The right panel shows the PDF weighed by mass, the differential mass distribution 
being the large plot, and the inset representing the cumulative distribution. It is clear that being able to trace a sizeable fraction of the 
baryonic mass (~ 10%), an instrument needs to go down to a surface brightness limit of a few CU. 



use for the detection algorithm, and thus a strong prediction 
of the simulation. We find a number density for the Lyman a 
bright sources of ri{Lya,z = 0.75) = 38 x 10"^ (Mpc/h)"^. 
Interestingly, the densities for bright sources in CIV and 
OVI are not very different : r]{CIV, z = 0.37) = 24.8 x 10"^ 
(Mpc/h)-^ and t){OVI, z = 1.1) = 17.3 x 10"^ (Mpc/h)-^ 
In these cases, however, the detection algorithms' param- 
eters do lead to a fraction of the faintest sources being 
dropped, i.e. it is not possible to find good combinations 
of the 2-fold threshold approach that isolates the point 
sources without introducing overlaps and/or losses. As 
those sources are too faint to be detected in any realistic 
scenario (see paragraph 6.2), we have not tried to fine-tune 
the detection algorithm in order to better isolate these 
objects. In any case, of much more importance, even for 
Lya compact sources, are the restrictions imposed onto us 
by the simulation's mass and spatial resolution. We have 
performed tests regarding conversion with earlier, lower 
resolution simulations, and found that we are 'complete' 
to a luminosity of log L(Lya)> 41.5 (at z=0.75; resp. 41.2 
for z=0.35), log L(OVI, z=l.l)^ 40.5, and log L(CIV, 
z=0.35)^40.0 (see below for calculation of the source 



threshold that could have been picked up for the more lenient 
cuts. 



luminosities). That is, increased resolution may result 
in producing a higher fraction of areas that had fallen 
underneath the threshold criteria used here in the lower 
resolution simulations, due to the strong non-linearity of 
the brightness as a function of density and temperature. 
Thus, for all transitions, the numbers quoted here are 
strictly lower limits. Furthermore, we remind the reader 
that we have extracted these sources in three different cubes 
for each transition that were not randomly chosen, but 
care was taken to select sightlines that do show large-scale 
structure. While this procedure resembles in some respect 
the target selection for a real observation, we nonetheless 
point out that cosmic variance is an additional factor in 
producing fluctuating number densities for which we have 
not accounted here. 



How big are the objects and what are their shapes ? 
The answers to these questions are not quite as straightfor- 
ward as one may naively assume. While the spatial extent 
perpendicular to the sightline in our cubes immediately 
translates to a size measurement, the coordinates along 
the sightline have the additional complication that the 
pixels we analyse have their peculiar velocities imprinted 
on them. As a flrst step towards estimating the sizes of 
the sources we have hence constructed small cutouts of 
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Figure 6. Same plot as Figure [S] but this time for OVI emission at z=l.l. In contrast to the Lya radiation, here most of the cube's 
volume is filled by gas in a regime (low temperatures, low densities and extremely low, if not zero, metallicity) that does emit very few, if 
any photons, hence the integral over the surface brightness range depicted here does not add up to unity. Note how much less important 
the cuts for self-shielding are regarding the distribution of the brightest voxels in the metal line, indicating that the main contribution 
to the flux in OVI (and CIV, see Figure (Tjl is emitted by gas with temperatures higher than the limits imposed by our self-shielding 
treatment. Note also how the some of the curves are exactly overlapping, as the self-shielding cuts do not affect any areas contributing 
to the surface brightness levels shown here. 
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Figure 7. Same as Figure[5l but for CIV at z=0.35. 



the cubes (centered on the coordinates found by Duchamp 
as central pixels) and summed them up along the full 
extent of the sources (given by their minimal and maximal 
estimate for the z coordinate in Duchamp). We then try 
to fit a simple 2-D Gaussian to the hght profile of such 
a collapsed image. In the majority of the cases the fitted 



FWHMs of the Gaussian (in x and y direction) are of a 
similar size, indicating that the underlying light profiles can 
be reasonably well described by the average of them for a 
crude size measure (Figure [TT] shows an example of such 
an object in Lyman a). For some objects, however, this 
blunt approach fails because their spatial (and often also 
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Figure 8. Similar to Figure[5l but now emphasising the redshift evolution for the Ly« emission from z=0.75 to z=0.35, and employing the 
two most extreme treatments for the self-shielding gas. The peak of the distribution shifts to lower surface brightnesses with decreasing 
redshift, mostly due to a weakening UV-background and dilution of the gas because of cosmic expansion. The bright, isolated areas 
show a slight increase in their brightness, but not as much as expected from a simple (1+z)'' scaling, rendering the evolution a more 
complicated function with redshift, density, temperature, and environment contributing to the complexity. 



velocity, see below) structure is more complicated. Figure 
[12] highlights a source that shows multiple bright cores, and 
even a hint for 'bridges' between them. 

Figure [9] shows the distribution of luminosities of the 
distinct sources detected by Duchamp in the Lya and metal 
line cubes, smoothed to FIREBall resolution, and according 
to the two most extreme treatments for the cells containing 
self-shielding or star- forming gas (CUTS 7 and 8). For 
Lyman a, the luminosities of the objects picked up by 
Duchamp is strongly dependent on the specific treatment 
of the observed cubes with regards to the self-shielding 
and star- forming g figure [5] clearly demonstrates. The 

various cuts introduced earlier lead to luminosities for the 
same objects straddling almost a factor of 100 in brightness 
in the extreme cases. Hence, for the Lya sources we can 
only bracket the luminosities assuming either the severe 
cuts (resulting in a median luminosity of log L [erg/s] = 
40.9 ) or the optimistic cuts (median log L [erg/s] = 41.8). 
The estimates for the OVI and CIV emission, on the other 
hand, are much more robust, because - as we have seen 
earlier - the cuts we introduce are not affecting pixels of 
the highest emissivity in these two transitions (median log 
L[erg/s] = 39.2 and 38.5 for OVI and CIV, respectively). 
Notice also, that the distributions for the bright objects in 
the metal lines are much broader than for Lyman a. This is 
not, as we will see in the next paragraph, due to a broader 
distribution in sizes of these objects, but a result of a larger 
spread in the peak brightness of the most dominant pixels. 
Overall, the objects observable in OVI and CIV are in 
general at least 1.5 orders of magnitude less luminous than 
the HI sources, but are spread fairly evenly over almost 3 
decades. Keep in mind, however, that for an assessment 



of their observability the much lower redshift for CIV and 
hence luminosity distance does help (yielding a factor of ~7 
in surface brightness). 

In Figure[TO]we show the relation between the luminosi- 
ties and the FWHMs (averaging the x and y components) of 
the compact objects - note that we fiagged out the obvious 
failures of the simple fitting. While the majority of the 
sources lie between about 50 and 100 kpc (proper), there 
are a few exceptions that can reach up to the enormous 
size of 300 kpc. The sharp low cutoff boundary (most easily 
visible for HI and OVI) represents the minimum FWHM 
of 2 pixels, i.e. is a result of the image resolutionlfl While 
there is a general trend for the more luminous objects in 
Lyman a to be also bigger, interestingly the reverse seems 
to hold for the CIV bright sources (and potentially also for 
OVI). Currently, we have no explanation for this anticorre- 
lation in size vs. luminosity for the metal line bright objects. 

When we instead collapse the 3D cutouts from above 
along the two spatial axes perpendicular to the line of 
sight, we can obtain some information on the velocity 
structure of the emitting gas. Figure [13] shows the profiles 
along the z direction for the a representative selection 
of bright sources in OVI along with their appearance in 
the spatial coordinates (green boxes) Here we plot the 
aggregated flux densities (integrated over the maxium 

^ Because of the lower redshift observable for CIV, the same pixel 
size in angular units span a smaller physical size, of course. 

The velocity profiles in CIV and Lya show very similar struc- 
tures, but we note that the interpretation of the Lyo structures 
is more complicated due to the variety of effects discussed earlier. 
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Figure 9. The luminosity distribution of sources extracted by Duchamp of cubes for three different transitions and redshifts, plotted 
differentially as the number density per (comoving) Mpc^h"^ in bins of log L = 0.25. The two different line styles for each transition 
are for the two extreme cuts for treating the self-shielding gas (CUT 8 being the most lenient, CUT 7 the most severe ; see text for 
details). Note that while the cuts have no effect on the total number of sources, a dramatic effect for the source luminosity in Lyo is 
apparent, whereas the distributions for the metal lines remain virtually unchanged. For Lya the cuts allow to bracket the most extreme 
cases. Notice how broad the distributions for the metal line transitions are in comparison to the rather sharply peaked hya luminosity 
distributions. This may partly be explained by compounding the complication of source 'incompleteness' due to the simulation resolution 
with the imperfections of the source extraction for the metal lines, whereas for Lyo we only need to to worry about the former (for 
details see text). 



source extent in x and y) versus the pixels' location in the 
cubes along z. If we assume the stretching in wavelength to 
be almost exclusively due to the peculiar velocities of the 
gas particles, we can translate the spectral extent into a 
velocity plot. Highlighted are four different scenarios into 
which we can group our objects. The zero velocity bin is 
defined arbitrarily to be the one with the highest flux, 
and we normalise by dividing through this peak flux. The 
top left panel shows the most common proflle : a single 
peak (representative of about 65% of the objects in Lya 
and OVI, but less than half in CIV - see table [3} . If we 
fit these single peak profiles with a Gaussian, the FWHM 
ranges between 150 and about 400 km/s, with the bulk of 
them between 200 and 300 km/sl^l The next largest group 
(25-30 % for all transitions) exhibits double profiles, which 
themselves can come in at least three different flavours : 
the majority has two peaks of almost equal height, and 

Keep in mind that one pixel of 0.25 A represents already about 
40 km/s. Hence, the lower limit is just a representation of our 
ability to resolve velocities. 



a trough in between that never reaches zero flux (right 
upper panel). Some of the objects, however, have one peak 
rather dominant over the other one (^ a factor of 2 in 
maximum flux). These can be further divided into sources 
where there is either no gap of zero flux in between them, 
or where there is a separation (but less than 200 km/s). 
Note that for most of these objects the spatial structure 
is still indicative of one single object, as in this example, 
but for a few sources it is clear that there are two distinct 
features that are responsible for the double peak in the 
veolcity proflle. A third category of objects has a multitude 
of different peaks, of varying heights with or without gaps 
of zero flux. While these represent less than 5% for the 
Lya and OVI emitting sources, they correspond to 10% of 
the objects for the CIV emission. While the velocity spread 
is in most cases for all of the transitions conflned to less 
than ±400 km/s from the zero velocity bin, there are a few, 
rare cases where the whole structure can span almost 1300 
km/s (see lower left and right panels). Those two examples 
again demonstrate that these complicated velocity proflles 
may be the result of many bright sources overlapping in 
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Figure 10. The size of the bright sources (measured as detailed in text via fitting a 2-D Gaussian to the images collapsed along z 
direction) vs. their luminosities. While there is a clear trend of more luminous Lya emission regions being more extended, the same does 
seem not to hold for OVI and CIV. The sharp cutoffs for the lower boundary of the sizes are a result of the minimum FWIfM allowed 
being 2 pixels, translating into different sizes at each redshift. Plotted are 3400+ sources in 3 different cubes for Lya, and 2600+ resp. 
600+ sources for OVI and CIV. All luminosities have been estimated using CUTS (the most lenient one), and utilise the results from 
the Duchamp source algorithm. 



the subcube extracted to derive the velocity profile (lower 
right), but need not be necessarily (lower left). While the 
former velocity profiles then can be readily explained by an 
overlap of a few single source profiles, the latter cases may 
be resulting either from multiple sources situated along 
the line of side, or be intrinsically more complex, possibly 
due to the imperfections of the extraction method (i.e. 
not taking into account spatial structures). Those latter 
cases, however, are only exhibited by a tiny minority of the 
sources (< 2% in all transitions). 



4.3 Comparison with Lya observations and 
origins of tlie bright Lya emission 

Our bright and compact sources, as displayed in Fig. 12 
and with properties summarised in table O reach Lya 
luminosities and space densities comparable to or larger 
than those obtained for the so-called Lya emitters (LAEs) 
revealed by a number of spectroscopic or narrow band 
imaging surveys at high redshifts (e.g. Ouchi et al. 2008 
and references therein). This comes as a surprise since these 
LAEs are identified as galaxies and their Lya emission 
essentially can be ascribed to star formation, a process not 



included in our emission models, but currently thought to 
be dominant in the Lya emission. How is that possible and 
why are the sources predicted by the simulations not seen 
in current surveys? 

Furlanetto et al. (20050 and Faucher-Giguere et al. 
(2010) have addressed these issues at 2; ~ 3. Here, we 
compare the luminosity distribution of the bright sources 
identified from our simulation with the only existing LAE 
survey at similarly low redshifts, the GALEX spectroscopic 
surveys at 0.20 < z < .35 (FUV) and 0.65 < z < 1.25 
(UV) ( Martin et al.l |2005| ). The GALEX surveys have 
found the LAEs to be less commo n and less luminous at 
the present time than at z > 3 |Deharveng et al. I I2OO8I : 
ICowie et al.ll2010l ). ICowie et al.l (|2011 ) have shown recently 
that most of this evolution occurs over the z = - 1 range. 
Since the luminosity distribution of the sources depends 
more stron g ly on the self shielding correction (see also 
lYang et al.l (|2006l )) than on the redshift, we limit ourself 
to the case of a single redshift, z=0.35, for which the 



Note however that Lya emission from star formation is in- 
cluded in their analysis. 
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Figure 11. Example of a compact bright Lya source at z=0.75 with a simple spatial profile (one single core ; center of image), as seen 
in a cube with a spatial resolution of 3 arcseconds . The image is a 40 by 40 pixel cutout, collapsed along the source's maximum extent 
in z direction. At this redshift, this translates to 800 (physical) kpc. The green contours delineate areas of +0.5 dex in brightness. Note 
how all three bright sources in this image are clearly separated, and are up to 4 orders of magnitude brighter than their surrounding 
areas. The unit of brightness here is LU rather than CU as we have collapsed along the spectral axis. 



luminosity functio n of LAEs is less uncertain than at 0.75 
l|Cowie et al.ll201ll ). At a Lya luminosity of about 3 xlO"*^ 
ergs , the lower limit for the LAE luminosity functions, 
the density of Duchamp sources (cut 7) is found to be 
23 times large r than the LAE density of 3x10"" Mpc~^ 
(AlogL=l)-^ (jCowie et al.ll2010t l^. At a Lya luminosity of 
- 10*^ ergs s"^ the excess factor decreases to 2 (for cut 7, 
but is as large as 100 for cut 8). 

There are however several reasons for missing objects 
like our simulated sources with the slitless GALEX spec- 
troscopic survey. First, a number of the sources from the 
simulation may not be associated to any continuum stellar 
light, whereas the extraction of GALEX spectra is based on 



Note that a comparison for cut 8 is less meaningful, since the 
curve is already declining at this luminosity value as a result of 
the method of analysis. 



the existence of a dispersed continuum. Second, in slitless 
spectroscopy line emission is smeared in the dispersion 
direction by the angular extent of the object. For a typical 
bright source size of 100 kpc (i.e. 20 arcsec at z = 0.35) 
the Lya FWHM would be increased from 7 A (spectral 
resolution) to about 20 A, resulting in a loss of contrast 
over the continuum and, at low signal to noise ratio, a 
non-identification as LAE or even a misidentification as 
AGN (if bright enough). Additionally, there is dilution in a 
direction perpendicular to dispersion, resulting in a loss of 
flux through the extraction window of 6 arcsec height. Even 
given all the uncertainties, we do not consider resonant 
scattering to be able to affect significantly the Lya photons 
from cooling radiation, in contrast to the situation with 
Lya photons resulting from star formation inside galaxies, 
where large optical paths through neutral gas may result 
in large quenching factors (if enc ountering dus t parti cles) 
as well as extended halos (e.g. ISteidel et al.l l|201ll) . In 
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Figure 12. Example of a compact bright hya source at z=0.75 with a more complex spatial profile. The image has the same dimensions 
as Figure rm Between the numerous sources in this image, there are 'bridges' connecting the bright sources that can reach up to 1/10 of 
the peak emissivity. 



Table 3. Properties of typical bright emission regions as detected with Duchamp, at the redshifts specified in the text (Ly« z=0.76, 
OVI z=l.l, CIV z=0.35). 



Transition 


Number density 
r?[#/(Mpc/h)3] 


Typical Maximum 
Surface brightness 


log L[erg/s] 
(median and range) 


Spatial extent 
in FWHM [kpc proper] 


Velocity structure 


Ly-a 


38 X 10-3 


> 3 X 10* CU 


41.8 {39..43.5) 
(lenient cuts) 

40.9 (39..42.2) 

(harsh cuts) 


majority : 50. .150 kpc 

extreme cases : 
> 200 kpc 


Single peak 60%, double peak 25% 
Multiple peaks : 5%, 

other : 10% 
Maximum extent for 
multiple peaks : 1300 km/s 


CIV 


25 X 10-3 


5 X 10^ CU 


38.5 (36..41) 
(cut independent) 


25.. 100 kpc 


Single peak 45%, double peak 30% 
Multiple peaks : 13%, 

other : 12% 
Maximum extent for 
multiple peaks : 800 km/s 


OVI 


17.3 X 10^3 


1.5 X lO'' CU 


39.2 (37.. 42) 


25.. 100 kpc 


Single peak 60%, double peak 30% 



(cut independent) 

Multiple peaks : 2%, other : 8% 
Maximum extent for 
multiple peaks : 1000 km/s 
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summary, the characteristics of the existing GALEX 
slitless spectroscopic survey may be able to explain the 
non-detection of extra Lya sources predicted by simulations. 

The fact that star-forming particles are excluded from 
our simulations, combined with the angular extension of our 
Lya Duchamp sources, suggest naturally an interpretation 
of these sources as CGM emission from Lya cooling radia- 
tion. This process has been repeated ly suggested to account 
for extended Lya emission (e.g . iHaiman et al.l (l200Ctl: 
iFardal et all l|2001^ : lDiikstra et al.1 (|2006l ): lDiikstra fc Loebl 
especially in the context of the so-called Lya blobs 
(e.g.lSteidel et al.l (|2000l ): lMatsuda et al.l (|2004l ): [Saito et al.l 
l|2006[ )'). This interpretation would be consistent with our 
predictions of a larger density of sources with respect to 
that of LAEs, since LAEs are themselves a small fr action of 
the st ar-forming galaxies (5 percent at z ~ 0.3, Co wie et al.l 
l|2010l ), 25 percent at 2 ~ 3. IShaplev et al.1 (j2003f ). for LAEs 
defined by an equivalent width W(Lya) > 20 A). 

The importance of Lya cooling radiation has been 
chall enged, however, by recent simulation s (jZheng et al.l 
I2OIOI ) and observations (|Steidel et aLllioiH ). showing that 
spatially extended Lya emission is a generic feature of 
high-redshift star-forming galaxies, resulting from Lya 
photons initially produced by star formation and reso- 
nantly scattered by neutral atoms in the circumgalactic 
medium. Because of surface brightness threshold effects, the 
extended emission would hence be observed individually 
only in a limited number of cases. This would lead to 
an underestimation of the total Lya flux by an average 
factor of 5 (Steidel et al. 2011). Such a factor, if applied 
to the luminosities of GALEX LAEs at z ~ 0.3, would 
be enough to make the resulting Lya luminosity function 
dominant over the luminosity distribution predicted from 
Lya cooling, at least for the cut 7 scenario (Fig. 11). In 
conclusion, this comparison validates our interpretation 
of the Lya bright sources extracted from our simulations 
as CGM emission, but leaves the contribution of cooling 
radiation to the total Lya emission observed within the 
uncertainties of simulati ons, and therein especially o f the 
self-shielding corrections l|Faucher-Giguere et al.|[2010l ). 



5 EXTENDED SOURCES/ WHIM 

This section gives a rough overview and some example plots 
for the structures that can be classified as 'filaments', i.e. 
here we take a look at the photons coming from regions 
outside the areas dubbed 'bright objects'. The main purpose 
is to be illustrative rather than quantitative and detailed, 
although some general features are summarised. 

The procedure to produce the graphs below is to use 
the cubes resampled to FIREBall resolution (as described in 
2.3), the specific cut values imposed are not important for 
these structures, the units of the data displayed here are in 
LU — #photons/(sec * cm^ * sr), as we integrate over the 
wavelength coordinate. The colour scaling of the images are 
chosen such that areas identified by the 3D-source finder 
Duchamp as 'bright sources' appear as white, compact 



'blobs'. The fiux outside these sources represents only a tiny 
fraction of the total flux (depending on the speciflc cut used 
less than 5% in the best case for Lya, and even less in the 
metal lines). In addition, much of the remaining fiux when 
flagging out the pixels identified as belonging to compact 
sources, is contained in areas very close to the volumes cut 
out (.i.e. < 50 kpc vicinity). Note that the contrast in flux 
density between the filaments and the bright blobs is up to 
4 orders of magnitude, whereas there is another factor of 
~50 in the surface brightness delineating the filamentary 
structures from the overall 'average' emission level. 

Fig. [2] shows a typical situation for Lya emission, 
as viewed by an instrument with roughly the specifics 
of FIREBall in angular and spectral resolution, but a 
much larger FOV (extended to 900 x 900 arcsec^). The 
filamentary structure connecting the bright sources has an 
extent in the x-y plane perpendicular to the line of sight of 
about 7200 kpc (physical), and extends over 6500 kpc in the 
z direction, assuming that velocity differences are not the 
main cause of the stretching over a wavelength range of 4.75 
A. Naturally, the boundaries of such structures are less well 
defined as the ones of the bright, isolated sources, so here 
we adopt the arbitrary definition whereby the outermost 
limit corresponds to the contour of the flux level 10 times 
higher than the median flux of all pixels. Then, the width 
of the long bridge, as measured perpendicular to the main 
direction between the bright sources, varies between 500 
and 750 kpc. This rough 10:1 ratio of length to width is 
quite representative of the fllamentary structures in Lya. 
An estimate by eye yields about 180 of such filaments in a 
cube of 25 x the FIREBall FOV, corresponding to a length 
per volume of 2.5 xlO^^ (Mpc/h)"'^, and a volume filling 
fraction of 0.1%. It is interesting to compare these estimates 
with the values of (Bond et al. 2010), who find and define 
filaments in their simulations via an algorithm that uses 
the eigenvectors of the Hessian matrix of the smoothed 
galaxy distribution. While their length measurements are 
of very similar size compared to ours (depending on the 
smoothing lengths they find l/V = 1.9 x 10"^ or 7.6 x 10"* 
(Mpc/h)^^), their filling factors (combining this length 
estimate with their mean filament widths) are of the order 
of 5%. This suggests that we are tracing similar structures 
(i.e. 'bridges' from and to massive and bright nodes of 
the Cosmic Web), explaining the similarity in overall 
length, but our focus on the HI emission outside of bound 
structures presumably results in much thinner sheets, as we 
are forced towards the densest parts of the gas. Following 
the definition above, the outermost regions of the filaments 
reach a surface brightness of about 0.1 CU, and even inside 
the filaments it usually remains below 5 CU, although very 
rarely bright spots (cf . the red- white areas of Figure I14|) can 
obtain up to 50 CUs. This is still a factor of at least 5 below 
the lowest brightness for the bright sources, but at least 3 
orders of magnitude brighter than the median flux density 
of the whole cube. Typical physical conditions of the gas 
encountered inside such filaments are roughly as follows : 
the density ranges from n= 5 to 50 x 10~^cm~'^, whereas 
the temperatures are between ~ 50, 000 < T < ~ 120, OOOK. 

Similar filaments can be seen in the metal line tran- 
sitions OVI and CIV. These still consist of 'bridges' 
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Figure 13. Examples of the 4 different classes of velocity profiles, here for the OVI emission. Note that the resolution of the representation 
chosen here (after resampling to the FIREBall specifics) is about 80 km/s per pixel. The coloured insets show the spatial images collapsed 
along the z axis as in Figure 1111 The velocity plots are derived by collapsing the area shown in green along the other two axis, and 
translating the differences in the z direction to velocities (for details see text). The upper left panel shows a single peak in the spatial 
image, but the z-direction yields a profile that can be fitted well by a single Gaussian. The panel on the upper right is an example of 
a similarly simply spatially structured source, but with a double peaked velocity profile. The lower left panel exhibits again a single 
peak in x-y, but this time a much more complex, multiply peaked velocity profile. The lower right panel, finally, has a complex spatial 
structure, in combination with a multiply peaked velocity profile. Note the slight change in scale for the velocity plot for this example. 
Both the CIV as well as the Lya compact sources show velocity profiles very similar to the ones depicted here, but the interpretation 
for the Lya emission is complicated by a variety of factors (see text for details). 



connecting two nodes (see Figures [15] and [16] for examples 
in OVI and CIV), however are much fainter than their 
Lya counterparts. In addition, the steeper decline of the 
emissivities in the metal lines with both temperature and 
especially density, leads to the remaining bright spots 
outside the compact sources to be more (OVI) respectively 
much more (CIV) concentrated around these sources than 
for Lya. By comparing the estimated gas emissivities in 
Lya, OVI, and CIV for this parameter space (see Figures [1] 
to [3| , and noting that they are all around the same value 
(assuming solar metallicity), we can conclude that the main 
reason for the metal line transition filaments to be so much 
fainter, is the gas metallicity, which tends to be slightly 
below 10""^ of the solar value, but can reach up to 10~^ in 
a few, isolated spots, usually near the bright sources, where 
it is reasonably to expect the metal enrichment of the IGM 
to take place first. 



6 OBSERVING STRATEGIES 



The idea of observing the emission from the IGM may be 
traced back to the early investigati ons of the natur e of th e 
UV (non-ionizing) backgro und (e.g. iDavidsen" et all ||1974 ): 
iParesce fc JakobsenI (|l98Cll ) and the idea that the resonance 
lines of HI and Hell from the IGM may contribute to this 
background. The role of clumping was alrea dy recognized 
as pr imordial for the possibility of detection (jParesce et al.l 
ll98d ). Since then, the association between the IGM emis- 
sion and the UV domain has come to the fo refront with the 
predictions of the WH IM and the CGM (|Cen fc OstrikeJ 
ll999l : [Dave et aLlliooH ). The development of this phase is 
expected to be maximum at low redshift and its tempera- 
ture implies observations only through Lya and high ion- 
ization species (OVI, CIV, etc) at UV rest-frame wave- 
lengths. Except for the Lya fiuorescence observations which 
make sense at high redshifts and can be attempted from 
the ground (|Rauch et al.l[200^ : ICantalupo et ah I I2OO5I ). the 
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Figure 14. Example of a typical filamentary structure seen in Lya emission at z=0.75. The size of the field in x and y is 900 x 900 
arcsec^, representing 6500 x 6500 kpc^ (physical) at this redshift. The slice seen here has a thickness of 5.0 A, corresponding to about 
6700 kpc in z direction. The bright, compact sources, which exhibit peak surface brightnesses of more than 4 orders of magnitude brighter 
than the filamentary structures connecting them, appear as white 'blobs' in the colour scheme chosen here. Note that the flux seen in the 
filaments is not spread uniformly, but strongly concentrated towards these bright sources, exhibiting a bright circum-galactic medium 
around these, before trailing off into 'bridges' forming the familiar sight of the Cosmic Web. Note that for Lya an additonal contribution 
to the surface brightness may arise from the resonant scattering of Lyman line photons (a.k.a as 'photon pumping', see text for details.), 
which is not included in this representation here. The physical conditions prevailing for the gas forming these large bridges are as follows 
: density n= 5 to 50 X lO^^'cm-^^ temperatures in the range of ~50,000 < T sS~120,000 K. 



WHIM and CGM-oriented observations have therefore to 
rely upon space or balloon experiments in the UV domain. 



6.1 Detectability of bright, compact sources 

What are the prospects of detecting the compact, bright 
sources identified via Duchamp as detailed in section 4.2. 
with current technology ? And what may be the best 
instrumental approach set out for this task ? Trying 
to gauge the detectability in a first approach, we have 
translated the Duchamp results into signal-to-noise-ratios 
(S/N) for each source, based upon the following simplifying, 
yet conservative assumptions. Note that certain aspects of 
the source detection already rely upon some fairly general 
instrumental characteristics (i.e. spatial and spectral reso- 
lutions). Furthermore, we assume to be dominated by the 
(cosmic) background as a noise source , and take this to 
be at a uniform level of about 500 CU dsrown et al.|[200ol : 
iMorrissev et al. l l2005l : lMurthv et al.|[201ff r For an effective 
telescope area of 230 cm^ (representative of a meter-class 
space-based UV-experiment) and an exposure time of 
1.0 X lO^s (representative of a typical deep exposure), we 
obtain the S/N ratios depicted in Figure[T7]for the three dif- 
ferent transitions at three different representative redshifts. 
The histograms have been computed from 3 emission-rich 
sightlines for each transition, totaling a combined 'observed' 



Field-of-View in each case of 3 xl5 x 15 arcmin'^. 

As can be easily seen, in this scenario we would be able 
to recover all of the Lya emitting sources, regardless of the 
specific cuts (solid and dotted black lines) or redshift (black 
or cyan curves) with high confidence (the vertical green line 
indicates a S/N of 5.0). Hence, even uncertainties in the 
precise source position and its extent, will not jeopardise 
detection at all - note that we have implicitly assumed 
to know both quantities perfectly in this approach for 
our simulations (degrading this precision by e.g. including 
random areas classified by Duchamp as not belonging 
to the compact sources only marginally affects the S/N 
for any sensible number of additional pixels, as we have 
tested). Thus, for Lya, the two main problems to solve are 
: a. guaranteeing a large enough FOV to allow for enough 
detected sources (e.g. to perform meaningful statistical 
analyses), and b. finding suitable tracers to guide target 
selection from observables established as good proxies 
to emission from the simulation. The latter point is, of 
course, only relevant for instruments with multi-object 
spectrographs, as sources with S/N >5 will be detectable 
by blind searches with integral-field spectrographs. 

The situation is somewhat different for the metal line 
transitions : both for the OVI and the CIV emitting regions, 
only a fraction of our simulated sources will be detected. 
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Figure 15. Example of a typical 'filamentary' structure in OVI at Ze,„ = l.l. The size of the field is again 900 x 900 arcsec^, translating 
to 7.4 X 7.4 Mpc^. The depth (5 A) represents 5.1 Mpc. Note the different scale in brightness as compared to Figure [Ml : the filaments 
connecting the bright sources are now about three orders of magnitude fainter. The physical conditions of the gas forming the filamentary 
structures are very similar to the ones seen in Ly« : the metallicity of such gas ranges from 10~^'^ up to 10~^'^, which is the main factor 
in those to be less bright. Note, however, that the flux is concentrated towards the bright sources even more than is the case for Lya. 



regardless of the specific cut used for the self-shielding, 
which is largely unimportant for the brightness of these 
systems, as detailed earlier. As the blue and red histograms 
in Fig. [nishow, about half of the CIV (at a redshift of 0.37) 
, and 2/3 of the OVI sources (at z=l.l) in our scenario may 
be missed. This fraction changes towards lower redshifts 
for OVI specifically, but note that the absolute number of 
sources per FOV then drops dramatically because of the 
enormously smaller area covered by the same opening angle 
(this effect can be seen e.g. in the two different redshifts 
for the Lya sources). Hence, an additional complication 
for targeting metal lines emission areas compared to the 
search for Lya sources is that one furthermore needs to 
be able to select tracers that guarentee to pick up the 
bright end of the distribution. Or be prepared to mostly 
rely upon stacking approaches because the majority of 
sources may individually lie below a detection limit, but 
collectively gather enough photons to statistically infer 
average quantities. 

Under the given assumptions the S/N scales with 



^Jtex-p X Af.ff. Hence for different combinations of exposure 
time and telescope size, the distribution simply needs to be 



shifted horizontally by the appropriate factor. Furthermore, 
the absolute source number is a linear function of the 
FOV, and hence the distribution can be shifted up- and 
downwards depending on this quantity, rendering this 
graphical representation a versatile tool to judge mission 
goals. 

A typical exposure time for a balloon-based experi- 
ment like e.g. FIREBall is about two orders of magnitude 
lower than indicated here - meaning that one can expect 
FIREBall to be successful in detecting Lya, and eventually 
a very small number of the brightest metal line emitting 
areas, depending on the size of the FOV. Note also, that 
in the case for Lya, we have specifically excluded photons 
stemming from star-forming regions, which may indeed be 
the major contribution to the fiux, even at larger distances 
from the galaxies themselves due to the resonant nature of 
the Lya transition (e.g. see Steidel et al. 2011 for examples 
at high z). While this additional contribution will facilitate 
detection of bright areas, it also complicates the analyses 
as the different components contributing to the total flux 
need to be separated. 
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Figure 16. Example of a typical 'filamentary' structure in CIV at z=0.35. The size of the field is again 900 x 900 arcsec^, equaling 
4.4 X 4.4 Mpc^, while the depth is 8.4 Mpc. The bright nodes are still connected by filamentary structures (with roughly the same 
physical conditions as in the OVI and Lya cases). 



Also note that we have thus far relied upon a fixed 
spectral coverage of only about 200 A, comparable to 
the capabilities of FIREBall. Extending this coverage by 
almost a factor of 10, as planned for certain space-based 
missions, will not only vastly increase the survey volume 
and hence the number of detectable sources - in this 
context it is worth reminding ourselves that we do not find 
a significant evolution in the Lya space density over the 
redshift intervals probed here, and thus the number counts 
should scale almost linearly with the survey volume - but 
also allow for the simultaneous observation of all three 
transitions, possible yielding important insights into the 
emitting material's physical state. 



6.2 Detection of filamentary structures 

As expected, a quick glance at the PDF for Lya emission 
(Figure [5] for Lya at z=0.75, but also valid for all other 
redshifts) confirms that by far the largest fraction of the 
cosmic volume lies below any realistic chance of a direct 
detection, as the peak of the distribution reaches mere levels 
of ~ 10~^ CU. Furthermore, a closer inspection of the flux 
remaining outside of the bright sources, which represent 
only 0.5 to 5% of the total flux to begin with, reveals that 



is is not homogeneously distributed, but rather shows a 
complex spatial structure in our simulation, with isolated 
spots reaching up to 50 CU over volumes as small as 5 voxels. 



We have conducted a few tests on a stacking approach. 
First, we have isolated by eye a few of the brighter fllamen- 
tary structures such as seen for example in Figure 1141 and 
added all their flux after removing the bright sources within 
them. Secondly, we have gone even one step further in our 
'naive' idealisation, and simply combined voxels outside the 
bright sources in descending order of brightness regardless 
of their actual position, until we maximised the S/N for the 
combined flux in that fashion (note that there is a maximum 
number of voxels after which adding new cells in fact leads 
to a diminished S/N due to the additional voxel being so 
dominated by noise rather than signal). Both methods, 
applied to a mock observation with an exposure time of 10® 
s and an effective telescope area of 230 cm^, lead to marginal 
detections of the aggregated voxels, irrelevant of redshift 
and the specific cuts for self-shielded gas. Again, recognising 
that these two scenarios are extremely optimistic or even 
overly idealistic - in a real observation e.g. there is no 
a-priori information where exactly the filaments themselves 
or their brightest part are located - it is obvious that dec- 
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Figure 17. The distribution of the signal-to-noise ratios expected for the objects found in the data cubes, as detailed in the text and 
shown in Figure[9l Here we are assuming an exposure time of 10® seconds, and an effective collecting area for an instrument of 230 cm^, 
typical for a proposed deep field observation in future space-based missions. The noise is assumed to be a combination of photon noise 
and the extra-galactic background, which we take as 500 CU. Note that the absolute number of objects is a function of the FOV and the 
spectral coverage. The latter has been chosen to be representative of FlREBall, and may be raised by a factor of 5 or more for future 
space-based missions. The former is in this case, a 3 X 15 X 15arcmin^ area of the sky. In addition to mass resolution effects (as is the 
main contributor for Lya) , the shape of the distributions for the metals is also influenced by the imperfections of the source extraction 
algorithm, especially for OVI towards the low S/N end. For details, see the text. 



tecting the faint, filamentary structure even just in Lyo0 
will remain a challenging endeavour with current technology. 

In reality, however, at least two aspects that may 
change this pessimistic outlook should be mentioned. First 
of all, as we have seen earlier, the mass resolution of 
the simulation may lead to an underestimate of the gas 
emission, specifically because it depends so non-linearly 
on temperature and densities. While this is most certainly 
the case for the brightest regions in the immediate vicinity 
of the densest and/or hottest areas, it seems less likely 
that large effects should be expected for the only slightly 
overdense and colder structures constituting the filaments, 
but we caution that even here clumping on scales smaller 
than resolved may occur. 



As we have seen in paragraphia there are indeed filamentary 
structures in the metal lines very similar to the ones observed 
for Lya, yet due to the metallicity being much lower than unity, 
these are so faint that we can safely focus on the - as we have 
shown - already difficult to detect emission from hydrogen alone. 



Secondly, it is worth comparing the emission calculated 
with our simulation, based upon photo- and coUisional 
ionization in an isotropic radiation field consisting purely of 
the average extra-galactic background, with estimates for 
at least one other process that we have identified already in 
2.4 as not being incorporated properly into our modeling, 
i.e. the photon-pumping of local non ionizing continuum 
photons from galaxies. In order to calculate the potential 
additional emission at the order-of-magnitude level, we 
have devised the following toy- model at z = 0.35: we have 
assumed all filaments to be dust-free uniform cylinders illu- 
minated by the UV continuum of young galaxies uniformly 
spread along their central lines. A 300 kpc proper FWHM 
diameter in HI has been taken, in rough agreement with the 
structures we have seen in paragraph 2.2. The cumulative 
comoving filament length per cubic megaparsec has been 
estimated by eye in one of our 900 x900 arcsec^ cones to 
be at ~0.006 Mpc~^. In the present coarse approximation, 
effects due to filament inclination on the line of sight would 
induce changes at the few tens percent level and have been 
neglected. Those assumptions result in an average column 
density along a filament radius of Ni^/ ~14.0 cm"'^, able to 
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scatter 0.4 A from a central continuum source when b = 
30 km/s, independently of the velocity distribution in the 
gas. The total luminosity scattered by HI in the filaments 
has then been taken equal to the total energy absorbed 
along all lines of sight, with the central source's luminosity 
chosen to produce the same luminosit y density near 1216 
A as measured by GALEX in the FUV (jSchiminovich et al.l 
[200^. We thus assume that most of the continuum near 
1216 A in a given portion of space is produced by flat con- 
tinuum galaxies distributed in the filaments, a reasonable 
assumption at moderate redshifts where the UV flux is 
dominated by galaxies of small mass. We find the HI Ly-a 
luminosity density "pumped" under those assumptions to 
be 4 times above that from our simulation. The average 
filament surface brightness is found to be of the order of 
10 LU, suggesting that areas with surface brightness in the 
hundred LU might not be uncommon given the luminosity 
spread of the illuminating galaxies in the real world, and 
the presence of cluste ring as already pointed out by e.g. 
iKoUmeier et all (|2010l '). This coarse computation brings the 
emission from such a filamentary structure into the realm 
of achievable sensitivities for a space-mission, provided the 
FOV is large enough to cover a significant area. Keeping in 
mind that velocity information may not be as crucial (or 
interesting) for this gas (cf. e.g. figure 1141 that shows the 
filaments not being spread over more than 100 to 200 km/s 
in velocity space), and also that high-spatial resolution is 
not mandatory, leaves us with the possibility of relaxing 
some of the instrument characteristics in order to maximise 
the S/N to be expected from the extended, faint emission 
regions. 

We caution, however, that while such an additional con- 
tribution to the hya radiation can boost the chances of ob- 
servability, it will tell us little about the physical state of 
the gas other than the distribution of neutral hydrogen and 
its velocity field (in conjunction with information about the 
non-ionising UV flux emanating from galaxies). The latter 
information though is relevant to the infall of cold gas. 



7 SUMMARY 

Exploiting a state-of-the-art large-scale AMR structure 
simulation, we examine the possibility of observing UV 
line emission at redshifts 0.3 ^ z ^ 1.2, focusing on the 
Lya and two metal line transitions (OVI and CIV). The 
emissivity estimates for the gas in the simulation utilises 
the spectral synthesis code CLOUDY, adopting options 
most suitable fo r the treatment o f the situation at hand 
(i.e. assuming a iHaardt fc Madaul l|200lj ) UV background 
with a zero percent escape fraction from galaxies as incident 
background at the simulation redshift ; turning off the 
line pumping for all transitions). Being interested in the 
emission from extra- and circumgalactic gas, we exclude 
all photons coming from cells in the simulation that are 
starforming. Furthermore, we are bracketing the emission 
of self-shielding gas by setting it either radically to zero 
emissivity, or by assuming coUisional ionisation equilibrium 
at the temperature and density provided by the simulation. 

Generally speaking, we encounter two types of emis- 



sion regions ; relatively compact, bright, isolated objects, 
and large, faint filamentary structures permeating all of 
the volume, connecting such dense knots. The former 
represent (presumably) circumgalactic material that is of 
high enough density and temperature in order to emit at 
surface brightnesses that can easily be reached with current 
instrumentation, while the latter will certainly remain 
challenging even for future satellite missions, as detailed 
below. By degrading the AMR simulation to the spatial 
and spectral resolution of FIREBall, an instrument with 
characteristics very typical for UV balloon observations, 
we can assess realistically the chances of observing the 
gas, and furthermore can describe in detail which phys- 
ical characteristics may be extracted from such observations. 

Applying Duchamp, a standard source extractor al- 
gorithm, to the observationally-minded datacubes, we are 
able to isolate the bright, compact objects for all three 
transitions at the three different redshifts dictated by the 
FIREBall spectral window (1990 - 2260 A, resuhing in 
0.64 <zem < 0.84 for Lya, 0.93 <Zem < 1.19 for OVI, 
and 0.28 <Zem < 0.46 for CIV). While extending over less 
than 0.1% of the cubes' volumes, these bright sources carry 
more than 95% of the flux, even for the least optimistic 
assumption for the treatment of the self-shielding gas. The 
number density for the Lyman a bright sources in our 
simulation is v{Lya,z = 0.8) = 38 x 10"^ (Mpc/h)"^, 
whereas the densities for bright sources in CIV and OVI 
are rj{CIV,z = 0.37) = 24.8 x 10"^ (Mpc/h)~^ and 
■n{OVI,z = 1.1) = 17.3 X 10^=* (Mpc/h)-='E!| While for the 
brightest areas in Lya the range depending on the specifics 
of the self-shielding treatment can be up to two orders of 
magnitude, the effect for the metal lines is minimal, as the 
main contribution to the light emitted stems in their cases 
from gas of temperatures well above any of the self-shielding 
cuts. Hence, our estimate for the median luminosity for the 
Lya sources ranges from log L [erg/s] = 40.9 to log L = 
41.8, while the median luminosities for the metal lines are 
robust against changes in the treatment of self-shielding 
; log L (OVI) = 39.2, and log L (CIV) = 38.5. Keep in 
mind, however, that while the luminosity distribution for 
Lya sources peaks sharply roughly at its median value, the 
distributions in OVI and CIV are much broader, without 
a clearly preferred peak. The overwhelming majority of all 
sources, Lya and metal lines, extend over roughly spher- 
ically symmetrical areas between 50 to 100 kpc (proper) 
FWHM (when fitting images collapsed along the velocity 
component/z direction of the cubes by 2D Gaussians), but 
in a few percent of the cases, the source can obtain sizes 
of up to 300 kpc. The morphologies of such large sources 
are often more complex, hinting at their being several, 
individually unresolved smaller ones in a group or cluster. 
Interpreting the sources' profile along the z axis of the cube 
as a velocity profile, we encounter a variety of different cases 
: on the order of 60% of the objects (again in all three tran- 
sitions) exhibit a single peak, with a median FWHM of 215 
km/s, and a range from 150 < a < 400 km/s. Next, there is 



^ The number of sources, however, is strongly dependent on the 
simulation resolution. See the Appendix for details regarding con- 
vergence. 
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a large group of sources (~25% of the cases) which show a 
double peaked profile, roughly split in half by ones that are 
symmetrical and others where one peak clearly dominates. 
Interestingly, for this group of objects with double-peaked 
velocity profiles their spatial profile looks in almost all cases 
still indicative of one single source, which is not true for the 
last group of objects with more complex, multiply peaked 
velocity structures. While some of the members of this 
category can have relatively simple spatial structure, many 
of them arc clearly examples of multiple, yet unresolved 
individual sources within the window over which we extract 
the velocity information. In general, the velocity spread 
over the whole source is less than ±400 km/s from the 
brightest pixel, but in rare cases can span almost 1300 km/s. 

In contrast to these bright, isolated sources with 
sharp 'edges', the morphology of the filamentary struc- 
tures cormecting the nodes of the Cosmic Web arc more 
complex with much less well defined boundaries. For all 
possible scenarios dealing with the self-shielding cuts, they 
collectively carry less than 5% of the flux in Lya, and 
even less in the metal lines. Adopting a straight-forward 
definition to delineate these from the surrounding medium 
(outermost contour has to have a surface brightness 10 
times larger than the median flux of the whole cube), 
we extract filaments that exhibit roughly a 10:1 ratio of 
width to length. The width, defined as length perpendicular 
to the longest segment connecting bright nodes, of these 
structures are between 500 and 750 kpc (physical). Within 
the definition above, the boundaries for the Lya filaments 
occur roughly at 0.1 CU, and remain inside those usually 
below 5 CU, although rare bright spots may obtain up to 
50 CU. The brightness of filaments in the metal transitions 
is about 2 orders of magnitude below that of Lya, largely 
an effect of the metallicity inside the filaments being of that 
order of magnitude below solar. Typical physical conditions 
of the gas encountered inside such filaments arc roughly 
as follows : density ranges from n= 5 to 50 x 10^''cm~'^, 
temperatures are between ~ 50, 000 ^ T ^ ~ 120, OOOK, 
and gas metallicity tends to be slightly below 10~^ of the 
solar value, but can reach up to 10^^ in a few, isolated 
spots, usually near the bright sources. The fiux distribution 
inside the filaments is not at all homogeneous, with about 
50% of it stemming from areas in the close vicinity of the 
bright sources (<~ 50 kpc). Given all these characteristics 
of the filamentary structures it is clear that their detection 
will remain a challenging task for the near future. 

Applying simplifying and conservative assumptions 
about source extraction as well as background and detec- 
tor noise, we conclude that the bright sources in Lya as- 
sociated with the CGM will easily be picked up with high 
confidence by instrumentation with current (or soon to be 
deployed) technology. For example, all of the simulated Lya 
sources at z=0.76, and a large fraction of the metal line 
transitions OVI and CIV exhibit S/N ratios in excess of 10 
when 'observed' by a meter-class space-based UV telescope 
with an effective area of 230 cm^ in a moderately long expo- 
sure of 10** seconds. While a ballon-based experiment with 
shorter duration observations will also be able to retrieve the 
majority of the Lya objects, the fraction of objects bright 
enough in the metal lines drops to the few percent level. A 



spatial resolution of 3-5 arcscconds suffices to minimally re- 
solve the bright CGM areas extending up to 200 kpc, and 
is currently easily obtainable technologically. Furthermore, 
a spectral resolution of ~0.5 A allows to extract the veloc- 
ity profiles of these objects, and hence obtain kinematic in- 
formation on the light-emitting gas distribution. The space 
density of the bright sources is high enough that even with 
a relatively short spectral coverage of ^ 200 A, already fea- 
sible with current detectors and dwarfed by possible space- 
based instruments that are poised to cover much more of the 
UV, surface densities can be achieved that guarantee suffi- 
cient numbers of objects in a survey field of a few arcmin^ , 
such that meaningful population analyses can be performed. 
Given these scenarios, the next question to be addressed has 
to be how to select good tracers in a possible surveying cam- 
paign. Note that we have purposefully neglected all stellar 
light emanating from the galaxies themselves, or even scat- 
tered Lya radiation, and hence our surface brightness esti- 
mates are strictly lower limits. In fact, distentangling such 
different components from the radiation we have simulated 
in our approach will become an important caveat for any 
survey. Hence, simply choosing UV-bright galaxies as tar- 
gets may not be the best way forward. Other source charac- 
teristics (e.g. metallicity, or environment) may serve better, 
and we are currently investigating this question. 
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APPENDIX : CONVERGENCE TESTS 

Especially towards areas of high density and/or temperar- 
ture, the ability to resolve small enough scales to capture 
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many of the relevant processes at least to some detail is cru- 
cial. As pointed out in the main text, e.g. the number of the 
bright sources is a strong function of the simulation (mass 
and spatial) resolution. A variety of tests can be performed 
in order to assess whether a simulation has converged onto 
stable values in a a wide range of parameter space. Because 
we are mainly interested here in the light emission of the gas 
outside startforming particles, we focus in this Appendix on 
two aspects regarding the bright, compact sources, and the 
surface brightness of all voxels outside these. 

Fig. [18] shows the luminosity distribution of bright 
sources in Lya isolated in exactly the same manner for 
two simulations with different resolution. Both simulations 
are examined at the same redshift (z=0.75), and over 
exactly the same volume. Whereas the black histogram 
shows the luminosities of detected sources in the original 
high resolution simulation, when applying the most lenient 
approach for handling the self-shielding, the red curve 
depicts sources found in a simulation that was run with 
half the resolution in each spatial direction, and hence 
with a mass resolution of a factor of 8 less. Clearly, the 
number of such bright sources is a strong function of the 
resolution : there are about 5 times as many such bright 
sources in the high resolution simulation. Ifowever, it is also 
obvious that almost all of the 'new' sources (as compared 
to the lower resolution output) are of lower luminosity, 
and their combined emission does not add substantial light 
to the total Lya radiation, which is dominated by the 
most luminous sources. Above an approximate limit of log 
L (Lya) ~ 41.6 both curves start to match. Hence, we 
conclude that the simulation has converged down to that 
limit of luminosity. In this sense, the number densities for 
the bright sources are strictly only lower limits. 

Outside of the bright sources, however, the situation is 
different. Fig. 1191 shows the surface brightness distributions 
of the same two simulation outputs mentioned above, after 
in each one the voxels identified to be belonging to bright 
compact sources (and a cushion of 5 pixels around these) 
have been removed. Those two histograms are remarkably 
similar, with only a very tiny minority of areas exhibiting 
a slightly brighter flux in the high-resolution cube. In that 
regard, we are able to assess that for the mildy over- to 
underdense regions singled out by cutting out the bright 
sources the simulation, not very surprisingly, has converged, 
rendering the predictions for those areas stable. 

Note that for both metal lines similar arguments hold. 

This paper has been typeset from a T^rjX/ I^Tp^jX file prepared 
by the author. 
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Figure 18. The luminosity distribution of the bright Lyo sources at z=0.75 extracted with the same source finding parameters in two 
cubes of the same volume with different spatial and mass resolutions. The black curve shows the results obtained in the simulation used 
throughout the paper, whereas the red histogram is derived from a simulation with a spatial resolution of a factor of 2 lower (in each 
dimension), and hence a mass resolution of a factor 8 worse. Note the strong dependence of the number of sources on the resolution 
: increasing the resolution in this fashion leads to about 5 times more sources in the same volume. The majority of the new sources, 
however, are much less luminous, and hence their contribution to the total light is small. The two curves start matching around log L ~ 
41.6, and hence we conclude that down to this luminosity the simulation has converged. 
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Figure 19. The surface brightness distributions of two Lya cubes at z=0.75 after removing the bright sources. The blaclt curve shows 
the results obtained in the simulation used throughout the paper, whereas the red histogram is derived from a simulation with a spatial 
and mass resolution of a factor 8 worse, as in Fig. 1181 Unlike for the bright sources, there is little difference in the two histograms, 
indicating - as expected - that for areas of moderate overdensity the simulation has reached convergence. There is a small number of 
pixels inside filamentary structures that exhibit surface brightnesses of about half a dex higher in the higher resolution cube, but those 
are so rare and spread out over wide volumes that they do not contribute significantly to the brightness of one specific filament. 



